Wormhole Effect in a Strong Topological Insulator 
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An infinitely thin solenoid carrying magnetic flux $ (a 'Dirac string') inserted into an ordinary 
band insulator has no significant effect on the spectrum of electrons. In a strong topological insulator, 
remarkably, such a solenoid carries protected gapless one-dimensional fermionic modes when $ = 
hc/2e. These modes are spin-filtered and represent a distinct bulk manifestation of the topologically 
non-trivial insulator. We establish this 'wormhole' effect by both general qualitative considerations 
and by numerical calculations within a minimal lattice model. We also discuss the possibility of 
experimental observation of a closely related effect in artificially engineered nanostructures. 

PACS numbers: 73.43.-f, 72.80.Sk, 73.20.-r 



Surface electrons in a strong topological insulator 
(STI) [THE] form a gapless helical liquid, protected by 
time reversal symmetry (T) through the topological in- 
variants that characterize the bulk band structure. When 
T is broken, which may be accomplished by coating the 
surface with a ferromagnetic film, the helical liquid trans- 
forms into an exotic insulating state characterized by a 
precisely quantized Hall conductivity 

with n integer. This result follows from the microscopic 
theory of the surface state [3J and also from the effective 
electromagnetic action describing the bulk of a topologi- 
cal insulator, which contains the axion term [9j [10] . Al- 
though it might not be possible to measure this 'frac- 
tional' quantum Hall effect in a transport experiment 
[TT] , Eq. |l]) is predicted to have observable physical con- 
sequences, such as the low-frequency Faraday rotation [9 
and the image magnetic monopole effect [12j . 

It is instructive to apply Laughlin's flux insertion ar- 
gument [T3] to the STI surface described by Eq. ([T]). 
This argument was devised to establish the fractional 
charge of quasiparticlcs in fractional quantum Hall liq- 
uids (FQHL) [14] and involves the adiabatic insertion of 
an infinitely thin solenoid carrying magnetic flux $(i) 
into the system, as illustrated in Fig. [T^,. As the flux 
is ramped up from to <f>o — hc/e, a circumferential 
electric field is generated in accord with Faraday's law 
V x E = — (l/c)(9B/9i). This induces a Hall current 
on the STI surface j = cr xy {E x z) which brings electric 
charge 

SQ = o~ xy — = ( n+ 2/ e ^ 

to the solenoid. Since the flux tube carrying a full flux 
quantum $o can be removed from the electronic Hamil- 
tonian by a gauge transformation, one concludes, as in 
FQHL, that an excitation with fractional charge |2]) must 
exist. This finding stands in contradiction to the well es- 
tablished microscopic theory of these surface states given 



by an odd number of massive Dirac Hamiltonians [3J . El- 
ementary excitations of a massive Dirac Hamiltonian are 
particle-hole pairs which are charge neutral. Yet, this 
same Dirac Hamiltonian exhibits Hall conductivity ([!]), 
which, through Laughlin's argument outlined above, im- 
plies fractionally charged quasiparticles. 

The resolution to this paradox comes from the realiza- 
tion that the quantum Hall state realized on the surface 
of a STI is inextricably linked to the bulk of the STI. 
Laughlin's argument fails because the flux tube inserted 
into the bulk of the STI is not inert. We demonstrate 
below that when $ = (s + l/2)<I>o, with s integer, the 
flux tube carries topologically protected gapless fermionic 
modes and forms a conducting quantum wire - a 'worm- 
hole' - along which the accumulated surface charge can 
escape to another surface of the sample. In the end, no 
net fractional charge is accumulated at the surface and 
Laughlin's argument instead predicts, indirectly, a new 
effect associated with a Dirac string in the bulk of a STI 
that we propose to call a 'wormhole effect'. 

In the rest of this Letter we establish the wormhole ef- 
fect, first by an analytical calculation using the universal 
properties of the surface states, and then by numerical 
calculations within a lattice model of a STI. We discuss 




FIG. 1: (Color online) a) Topological insulator coated with a 
ferromagnetic (FM) film. The flux tube employed in Laugh- 
lin's argument and the induced electric field are indicated, b) 
Flux tube threading a cylindrical hole in a STI. Arrows illus- 
trate the helical spin state for upward moving electrons (for 
down-movers the arrows are reversed). 
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its physical properties, significance and the possibility of 
experimental observation. 

We begin by considering a bulk STI with a cylindrical 
hole of radius R threaded by magnetic flux $ = ?7$o with 
< 77 < 1 as illustrated in Fig. [TJd. By solving the Dirac 
equation for the surface electrons we show that a gapless 
state exists when r\ = | and persists in the limit R — > 0. 
According to Ref. [TS] electron states on a curved surface 
of a STI, characterized by a normal unit vector n, are 
described by a Dirac Hamiltonian of the form 

H = ^v[HV- n + n- (p x a) + (p x a) ■ n] (3) 

where v is the Dirac velocity, p = ~ih\J is the momen- 
tum operator and cr = (ci, 02, 03) is the vector of Pauli 
spin matrices. The magnetic flux is included by replac- 
ing p with 7r = p — (e/c)A, where A — i]$o(z x r)/2irr 2 
is the vector potential. For a cylindrical inner surface 
n = — (cos tp, sin if, 0), the Hamiltonian Q becomes, in 
cylindrical coordinates and taking v = H = 1, 



H k = - 



1 

2R 



-ike~ tv 



(4) 



We assumed a plane- wave solution e tkz along the cylinder 
axis and replaced — id z — > k. 

The eigenstates of Hk are of the form 



*fc(¥>) 



(5) 



with I integer. The spinor VE^ = (/&, <?fc) T is an eigenstate 
of T-Lki = 02 & — 03 (i + I — rf) /R with an energy eigenvalue 



E u = ±vh 




(6) 



For a generic strength of the magnetic flux the spectrum 
of electrons along the cylindrical surface shows a gap 



A 



2vh 
~R 



(7) 



When rj — |, i.e. at half flux quantum, the I = mode 
becomes gapless, E^q 1 ^ 2 = ±vh\k\, independent of the 
hole radius R. This is the wormhole effect introduced 
above. Physically, the necessity of the flux for the gapless 
state to occur stems from the Berry's phase n acquired 
by electron spins in the helical state depicted in Fig. [T]d. 
The gapless state occurs at half flux quantum when the 
Aharonov-Bohm phase exactly cancels the spin Berry's 
phase. 

We observe that the system remains T-invariant in 
the presence of a half flux quantum threading the hole. 
Therefore, the gapless state is topologically protected 
against any weak perturbation that respects T and does 




FIG. 2: (Color online) Band structure of 36 x 18 STI sample 
(periodic boundary conditions), infinite along the z-direction 
with two flux tubes carrying flux t;"I>o and — 77<I>o inserted, a) 
77 = 0, b) 77 = 1/2. Panel c) shows dependence of the spectral 
gap on the flux for thin flux tubes threading a lattice plaquette 
(solid line), and thicker flux tube threading a 3 x 3 rectangular 
hole (dashed line), d) Wavefunction amplitude for a low- lying 
state extended along the flux tubes. Parameters used here are 
t = 0.2, e = 0.8. 



not close the bulk gap. Specifically, it should be ro- 
bust against weak non-magnetic disorder as well as any 
smooth deformation of the hole. Our numerical simula- 
tions, presented below, provide support for this topolog- 
ical protection. 

We now study the wormhole effect using a concrete 
lattice model of a topological insulator which we solve 
by exact numerical diagonalization. In order to keep 
the computational difficulties at a minimum we consider 
a simple model on a cubic lattice discussed previously 
[21 [TH]. This minimal model has two electron orbitals 
per lattice site, denoted c and d, and is defined by the 
momentum space Hamiltonian H — ^ k 4'j c 'Hk , I , k with 
*k = (ckt ; c k ;,dkt,rfH) T , 



-2A T z cr M sin k p 



T x m k , 



(8) 



and mk = e — 2tY] cos fc^. Here and are Pauli 
matrices in orbital and spin space, respectively, with 
fj, = x,y, z. The system defined by H is invariant un- 
der time-reversal and spatial inversion. The spectrum of 
excitations has two doubly degenerate bands, 



£k = ±\/4A 2 (sin^ k x 



sin 2 k y 



sin 2 k z ) 



mi 



(9) 



At half filling, depending on the values of the parameters 
A, t, e, the system can be a trivial insulator, as well as 
a strong and weak topological insulator (WTI) [HI fTB] . 
Below, unless stated otherwise, we work with parameters 
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2t < e < 6t, corresponding to a STI phase characterized 
by the Z2 invariant (1;000). All energies are expressed 
in units of A which we take equal to 1. 

To look for gapless propagating modes along a flux 
tube we first consider a sample infinite in the z-direction 
with a rectangular base containing 2L x L sites. Two 
straight flux tubes carrying fluxes rj&Q and — along 
the z-direction are positioned a distance L apart on the x 
axis. Since the total flux threading the system is zero for 
this arrangement we may use periodic boundary condi- 
tions along x and y and thus eliminate gapless modes that 
would otherwise reside on surfaces. Results for L = 18 
are displayed in Fig. [2] Without the flux (rj = 0) the 
system shows a spectral gap. For 77 > subgap states 
appear near k = 0. When 77 = \ a gapless mode exists 
along each of the flux tubes. 

In addition to thin flux tubes that thread an elemen- 
tary plaquette we also considered flux threading a larger 
rectangular hole in the sample. Results for this case are 
similar; a gapless mode appears when r\ = | and the 
dependence of the gap on 77 (Fig. [2J;) now more closely 
resembles that given in Eq. ^ . In general gapless modes 
persist for any size and shape of the hole, as long as it is 
threaded by a half flux quantum. 

We performed similar calculations for other topologi- 
cal phases occurring in the same model. In a (1; 111) STI 
that occurs when — 6t < e < — 2t a single gapless mode 
(per flux tube) exists, now located near k = tt. In WTI 
phases an even number of gapless modes per flux tube 
appear. For a straight flux tube along direction n we 
find two gapless modes (one at k = and one at k = 7r) 
when n • v ^ and zero otherwise. Here v = {v\Viv-s). 
Finally, we have verified that in a trivial (0; 000) insula- 
tor, that occurs for |e| > 6t, no gapless modes appear for 
any direction or strength of flux tube. 

Using our model with a single flux tube in a geometry 
with open boundary conditions it is possible to visualize 
the flow of charge at the intermediate steps of Laughlin's 
flux insertion argument. To this end we consider a cube 
of size L 3 and supplement the Hamiltonian Q with a 
surface magnetization term 

h s = -ft s E • (*5**i) • ( 10 ) 

j£surf 

Here ij represents the unit vector pointing outward from 
the origin located at the cube's center and fis is the sur- 
face magnetization strength. Hs breaks T at the sample 
surface and a gap of size ~ 2|f2s| opens up in the spec- 
trum of the surface states. Figure [3^, shows the evolution 
of charge SQ accumulated near the intersection of the 
flux tube with the magnetized surface as a function of 
i] = &/&o- For small X] we observe SQ = \er\, consis- 
tent with the fractional Hall conductivity a xy = e 2 /2h 
expected on the basis of Eq. |lj. At rj = | a charge 
e/2 travels along the flux tube and combines with the 



a) 




FIG. 3: (Color online) a) Charge SQ (in units of e) induced 
in the lower half of the sample by a flux tube carrying flux 
$. Insets show the charge density at indicated values of 
flux, b) Charge SQ in the sphere of radius ro centered at 
the monopole as a function of monopole strength. We have 
chosen ro = 9.097 so that the sphere is large enough to enclose 
the expected charge e/2 for the unit monopole projecting flux 
"I>o- In both panels a cube-shaped sample with 20 3 lattice sites 
and parameters t = 1, e = 4 is used, Qs = 1,0 in a) and b) 
respectively. SQ is defined relative to the 77 = situation. 

negative charge that has built up on the opposite sur- 
face. For ij > \ the charge SQ grows again at the rate 
controlled by a xy until it reaches SQ = at rj = 1. As 
already mentioned above, a Dirac string carrying a full 
flux quantum <j> can be removed by a gauge transforma- 
tion and the above evolution is thus consistent with the 
expectation that this weakly interacting system returns 
to the original configuration at the end of a full cycle. 

Figure displays a modified arrangement with the 
flux tube terminated by a magnetic monopole located at 
the center of the sample. This furnishes a realization of 
the Witten effect [T7] in a STI [TB]. As a function of 
increasing rj, the charge first accumulates at the intersec- 
tion of the flux tube and the surface. At r\ — h a charge 
e/2 travels along the wormhole to the monopole, the cor- 
responding charge density clearly visible in the inset to 
Fig. [3)3. At r\ = 1 the flux tube becomes invisible but the 
e/2 charge remains bound to the monopole as expected 
on the basis of general arguments [16l [IT] . 

We performed similar calculations for flux tubes of var- 
ious shapes and in the presence of weak non-magnetic 
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disorder. In all cases we found low-energy modes associ- 
ated with the r\ = j flux tube confirming the topological 
robustness of the wormhole effect. 

We now address the possibility of experimental detec- 
tion of the wormhole effect predicted in this Letter. In 
a real physical system it is not possible to confine mag- 
netic flux to an area of size comparable to the crystal 
lattice spacing as would be necessary to probe the worm- 
hole effect in its pure form. However, it should be pos- 
sible to observe a closely related effect in a nanoscale 
hole fabricated in a STI crystal with a uniform magnetic 
field applied parallel to its axis, Fig. [TJd. Sweeping the 
magnetic field strength will result in a periodic varia- 
tion of the conductance along the hole with minima at 
n$ an d maxima at (n + l/2)<& as the excitation spec- 
trum oscillates between insulating and metallic. Such 
variations should be observable experimentally if certain 
conditions are met. First, the hole radius R must be 
sufficiently large so that several oscillations can be ob- 
served in the available range of the laboratory field B. 
This gives 

R <; (NQo/nB) 1 /' 2 for N oscillations. Second, 
R must be sufficiently small so that the maximum spec- 
tral gap Eq. Q is large compared to fc^T, or otherwise 
the oscillations in the conductance will be washed out by 
thermal broadening. This gives R < hv/V2k B T. Taking 
typical values B = 10T, N = 10, v = 5 x I0 5 m/s and 
T = IK yields 36 nm < R < 450 nm. Thus, the exper- 
imental challenge would lie in fabricating a sub-micron 
size hole (or an array of holes) in a STI crystal or a thick 
film and measuring the conductance along the holes. 

Oscillations with period $0 have been observed in re- 
cent conductance measurements on Bi2Se3 single-crystal 
nanoribbons (cross sections 6 — 10 x 10 _15 m 2 , consistent 
with the above bounds on R) in longitudinal magnetic 
field [23]. In these, the same effect as discussed above 
should occur for the topologically protected states on the 
outer surfaces of the nanoribbon. However, the observed 
positions of minima and maxima were opposite to those 
predicted by our theory, suggesting that conductance in 
these experiments is dominated by some competing ef- 
fect. The oscillations reported in Ref. [23] clearly deserve 
a detailed theoretical study. 

The wormhole effect introduced here is fundamen- 
tally different from the gapless modes predicted to exist 
along the core of a crystal dislocation in topological in- 
sulators [18]. These latter modes depend solely on the 
weak invariants {y-\y?y-&) whereas the wormhole effect de- 
pends on the more robust strong invariant Uq. In this 
sense the wormhole effect is inherently three-dimensional 
while the gapless modes associated with a dislocation are 
more closely related to the zero-modes in 2-dimensional 
topological (spin-Hall) insulators with solitonic defects 
P3H [2D] . Mathematically, the wormhole effect is related 
to the protected one-dimensional modes predicted to ex- 
ist along vortex lines in the order parameter characteriz- 
ing a topological Mott insulator [21]. Unfortunately, no 



topological Mott insulators are known to exist at present, 
although according to Ref. [22] the pyrochlore compounds 
A2IX2O7 (A — Pr,Eu) may exhibit this behavior. 

The wormhole effect studied in this Letter represents a 
distinct bulk manifestation of the unusual electron prop- 
erties in a strong topological insulator. Its existence re- 
solves a conceptual dichotomy that arises when Laugh- 
lin's argument is applied to the magnetized STI surface 
and exemplifies a unique bulk-surface correspondence in- 
herent to STIs. A closely related counterpart of the 
wormhole effect should be observable in artificially en- 
gineered nanostructures fabricated from available STIs. 
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Note added — When this work was close to completion 
we became aware of a preprint [15] which has identified, 
in a different context, the gapless modes on a surface of 
a STI cylinder threaded by flux <&o/2. 
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